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Abstract 

We report a first-principles projector augmented wave (PAW) study on Nao.5Co02. With the 
sodium ion ordered insulating phase being identified in experiments, pure density functional cal- 
culations fail to predict an insulating ground state, which indicates that Na ordering alone can not 
produce accompanying Co charge ordering, if additional correlation is not properly considered. At 
this level of theory, the most stable phase presents ferromagnetic ordering within the C0O2 layer 
and antiferromagnetic coupling between these layers. When the on-site Coulomb interaction for Co 
3(i orbitals is included by an additional Hubbard parameter U , charge ordered insulating ground 
state can be obtained. The effect of on-site interaction magnitude on electronic structure is stud- 
ied. At a moderate value of U (4.0 eV for example), the ground state is antiferromagnetic, with a 
Co^"^ magnetic moment about 1.0 and a magnetic energy of 0.12 eV/Co. The rehybridization 
process is also studied in the DFT-I-U point of view. 

PACS numbers: 71.28.+d, 71.27.+a, 75.25.+Z 



I. INTRODUCTION 



Naa;Co02 was originally studied as a kind of rechargeable battery material^ and ther- 
moelectric material2*2ii. More recently, the discovery of the superconductivity in its hy- 
drated x=0.35 compound^, where the effect of hydration is found to be limited to lattice 
expansion^"^, makes it receive a renewed interest. Nai.Co02 has a crystal structure consisting 
of 2D triangular lattice Co sheets, octahedrally coordinated with O above and below the 
Co planes, and layers of Na ions sandwiched between the C0O2 sheets. The phase diagram 
of Na2.Co02 has been determined by changing the Na content x using a series of chemical 
reactionsi^. Their electronic and magnetic properties are found to be strongly dependent 
on the number of charge carriers introduced by deviations from stoichiometry in the Na 
sublattice. 

Theoretical studies on the electronic structure of Na2^Co02 are very active. The early 
work by Singb^iiSifocused on the magnetic properties. He predicted a weak instabilities of 
itinerant ferromagnetic character in range x = 0.3 to x = 0.7, with competing but weaker 
itinerant antiferromagnetic solution. He used the virtual crystal approximation, instead of 
supercell model, to describe the partially occupation of Na sites. As a result, the possibility 
of sodium ion and Co^^/Co^"*" charge ordering is not allowed. Using pseudopotential plane 
wave method, Ni et al~ have optimized the geometries for Nai,.Co02 {x = 0.25, 0.5, 0.75 
and X = 1) with a 2 x 2 x 1 supercell. They found that the itinerant ferromagnetism is 
strongly suppressed by the local distortions of the oxygen around the cobalt. With optimized 
geometries, they identified a phase transition from wide-band ferromagnetic to narrow band 
paramagnetic metals with x increasing. By carefully examining the doped electron density, 
Marianetti et al.^ found that the rigid band model is not suitable for this Na doped system. A 
rehybridization driven by a competition between the on-site Coulomb interaction and the Cg- 
oxygen hybridization is identified. They have also used a modified Hubbard model to study 
this kind of rehybridization. Within the framework of density functional theory (DFT), the 
on-site Coulomb correlation effects can be dealt with DFT+U methocii^. Using this method, 
Zou et alM have studied the electronic structure of Naa;Co02 with varying c/a ratio. They 
did not find critical change of the band structure near Ep, except a notable decrease of the 
DOS at Ep. To consider the relations between correlation and charge ordering, Pickett and 
coauthorsi^ii^ have applied the DFT+U method to a \/3 x -\/3 x | supercell of Naa;Co02 
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at X = 1/3 and x = 2/3. They found that the parameter U is critical to charge and 
spin orderings (for small supercell, the so called charge disproportionation always leads to 
charge ordering). They suggested that the phase diagram of Na2;Co02 is characterized by 
a crossover from effective single-band character with U ^ W for x > 0.5 into a three-band 
regime for x < 0.5, where U ^ W and correlation effects are substantially reduced. 

Nao.5Co02 is of special interests among the NaxCo02 family. In the recent experimental 
phase diagram of Naa;Co02, a novel sodium ion ordered insulating state at a; = 0.5 was 
found^, with the most strongly developed superstructurei^^ii. This insulating state separates 
two kinds of metallic phase, i.e. Curie- Weiss metal for x > 0.5 and paramagnetic metal 
for X < 0.5. It is suggested that the occurence of the insulating state indicates a strong 
interaction between the ions and holes even though they occupy separate layers, and a charge 
ordering is thought to accompany the sodium ion ordering. Therefore it is interesting to 
compare the electronic structure based on this experimental sodium ion ordered structure 
with previous theoretical results to identify the effect of Na ordering. For example, it will 
be interesting to see if the Na ordering will produce a Co charge ordering at DFT level. 
On the other hand, the correlation effects of Nao.5Co02 is considered to be a betweenness 
comparing to small values for x > 0.5 and relatively large values for x < 0.5^^'^^, which 
demands a study at the DFT+U level to identify the effect of correlation. Unfortunately, 
there is no theoretical work in the literature for this important x = 0.5 case with both 
Na ordering and correlation effects being properly considered. In this paper, we report the 
electronic structure of Nao.5Co02 focusing on the possibility of charge ordering and based 
on the experimental superstructure at both DFT and DFT-I-U level. 

II. COMPUTATIONAL METHOD 

The calculations in this work were performed with the Vienna Ab initio Simulation Pack- 
age (VASP)^^*^, which is a first-principles plane-wave code, treating the exchange and corre- 
lation in the DFT scheme. In this study, the Perdew-Wang functional forn>^ of generalized 
gradient approximation (GGA) was used. For spin-polarized calculations, the spin interpo- 
lation of Vosko et al."^^ was also adopted. The projector augmented wave (PAW)^ method in 
its implementation of Kresse and Joubert^'^ was used to describe the electron-ion interaction. 
The Kohn-Sham equations were solved via a Davidson-block iteration scheme^^. 
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A -y/S X 2 X 1 supercell was introduced to describe the experimental sodium ion order- 
ing. Brillouin zone (BZ) integrations were performed on a well converged grid of 8 x 8 x 4 
Monkhorst-Pack24 special points. The total energy and density of states (DOS) were calcu- 
lated using the linear tetrahedron method with Blochl corrections^^. The plane wave kinetic 
energy cutoff was fixed to 500 eV. Dudarev et al. 's approach^^ was adopted for DFT+U 
calculations. Since this DFT+U functional depends only on the difference of Hubbard pa- 
rameter U and screened exchange parameter J, J was fixed to 1 eV during all DFT+U 
calculations. The Hubbard on-site Coulomb interaction is applied to Co 3(i orbitals only. 
Geometry relaxation is performed at a spin unpolarized GGA level. 

III. RESULTS AND DISCUSSION 
A. Geometrical properties 

A X 2 X 1 superstructure of Nao.5Co02 has been found by powder neutron diffraction 
experimenlfi^. This superstructure has the Pnmm space group symmetry, with a cell for- 
mula: Na4Co80i6. In this experimental structure model, the Nal and Na2 sites (26 and 2d 
Wyckoff sites for original PG^/mmc space group, and 2b and 2a sites for this Pnmm space 
group respectively) are equally occupied, and the ordered Na ions form one-dimensional 
zigzag chains. Two types of Co ions (4/ and 4d), which differ subtly in their coordination 
by oxygen, are also in chains. The existence of two kinds of Co sites allows the emer- 
gence of two inequivalent Co ions, namely Co^"*" and Co''"'", in the process of self-consistency. 
Therefore this structure model can be used to study the possibility of charge ordering. 

The geometry optimization, starting from the experimental structure, reaches its con- 
vergence until all forces vanished within 0.01 eV/A. The optimized geometry is shown in 
Fig. ^ During the optimization, the cell parameters are fixed to the experimental values 
(a = 4.876 A, b=5.630A, and c = 11.063 A), which results in only a neglectable stress less 
than 20 kB. As listed in Table HI relaxations of atomic positions from the experimental ones 
are also found to be very small. The largest relaxation comes from the first two Na atoms 
along X direction (0.06 A), typical relaxations are about 0.01~0.02A. The calculated two 
kinds of Co-0 distances (1.91 A and 1.88 A) are very similar to experimental values (1.90 
A and 1.86 A). In an ideal structure model of Na2,.Co02, Co and O may form layers of 
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edge-shared CoOe octahedra in a triangular lattice, but distortion of the octahedra with 
the variation of the O hight zq from its ideal value is always exists in real system. Our 
optimized zq varies from 0.0845 c to 0.0892 c for different O sites, resulting in the Co-O-Co 
angles varying from 96.2° to 97.0°. Comparing to the 90° angle in the undistorted octahe- 
dra, the optimized structure shows a considerable distortion of the oxygen octahedra. For 
2x2x1 supercell, we get a slightly smaller distortion, but Ni et al.^^ get a significantly 
larger one. We notice that they have used a relatively small kinetic energy cutoff in their 
study. 

B. GGA electronic structure 

We can clearly identify three well separated manifolds in the paramagnetic GGA band 
structure and DOS of Nao.5Co02, as shown in Fig. |21 The lowest in energy is the O 2p 
bands, with a bandwidth about 5.5 eV (only part of them are shown in the figure). The 
other two manifolds come from the Co 3d orbitals, which are splitted into a lower lying 

manifold and an upper lying Cg manifold, separated by approximately 2.5 eV according 
to the oxygen octahedral crystal field. Although these three manifolds are well separated 
in energy, the hybridizations of O 2p states with both the and Cg states are notable. 
The Fermi energy locates at the upper edge of the t2g manifold, indicating a metallic GGA 
ground state. 

The low energy properties of Nao.5Co02 will be determined by the t2g manifold, which has 
a bandwidth of 1.6 eV. The trigonal symmetry of the Co sites further splits the t2g states into 
one with ag symmetry and a doubly degenerate e'g pair. In Fig. |2l we also plot the ag (Co 
(1^2) partial DOS. As found by Lee et al}^, the ag states have a bandwidth almost identical 
to that of the whole t2g manifold, which indicates that the widely used single-band model 
may not be sufficient enough to describe the GGA electronic structure of this material. 

Fermi surface often plays an important role in the electronic structure of materials. Our 
GGA Fermi surfaces for the five bands across Fermi level are shown in Fig. 01 which are 
obtained by a 32 x 28 x 14 BZ sampling and viewed by XCrySDen package^l. LDA/ GGA band 
structure is expected to give valuable and accurate information about the Fermi surface even 
in system with strong correlations^. In our supercell structure, there will be band folding 
from the first BZ of the original hexagonal lattice (black hexagon in the Fig. ^^), and the 
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Fermi surface thus becomes very complex. We notice that, on the other hand, the band 
structure of the superstructure geometry is far from simply band folding, for example, the 
Fermi surface of the fourth band (band 99) shows a prominent 3D character, which has 
not been reported in previous virtual crystal study^^. According to the complexity of the 
present Fermi surfaces, nesting effect for Nao.5Co02 may not be very strong, and we thus 
don't expect to find corresponding spin density wave (SDW) state. 

Spin-polarized calculations are carried out in order to address the possibility of magnetic 
ordering. When the ordering is constrained to be ferromagnetic (FM), a stable low spin 
solution is obtained, with magnetic moments of 0.4 and 0.45 fiB for Col and Co2 respectively. 
But the corresponding energy gain from ferromagnetic instability is only 26 meV/Co. The 
moment values are somewhat smaller than that (0.5 /i^) reported by Singh^ under virtual 
crystal approximation, which may due to the difference between their structural model and 
ours. In our structure, the triangular symmetry is slightly broken. 

Antiferromagnetic (AFM) ordering is frustrated on the 2D triangular lattice. But if the 
magnetic moments on a sublattice vanish, the left other sublattice may geometrically permit 
the presence of AFM ordering. Here we get an AFM model (AFMl) by fixing the magnetic 
moments at Col sites to be zero. As shown in Fig. ^ the Co2 sites form a ID chain in the 
x-y plane, and it is easy to establish AFM ordering within this chain. In a unit cell, the two 
AFM chains are set to be out of phase. The AFMl model results in a magnetic moments 
of Co2 ions about 0.39 fiB, which is much larger the result (0.21 fis) for the Singh's AFM 
model^. We notice that in that model each Co ion has also two nearest neighbors of like spin 
in addition to four of opposite spin. Another AFM model (AFM2) is constructed by alter 
the direction of magnetic moments of neighboring FM Co layers. Although the interlayer 
exchange interaction is expected to be small, we get lower energy for AFM2 state than that 
of FM state. 

Until now, we get a metallic ground state with intra-layer FM ordering and inter-layer 
AFM ordering. Comparing with the results of Singb^i, we may conclude that Na ordering 
alone will not qualitatively affect the GGA electronic structure to generate the band gap 
and charge ordering observed in experiments. For the AFMl state, where charge ordering 
is enforced, the magnetic energy (7 meV/Co) is much smaller than those of AFM2 and 
FM states, which indicates the enforced charge ordering under GGA level is unfavorable in 
energy. This result is consistent with the calculations for Nai/3Co02, where attempts using 
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LDA to obtain self-consistent AFM spin ordering always converge to the FM or nonmagnetic 
solutionis. This discrepancy between theory and experimental observation naturally leads 
to a further consideration of correlation effect. 

C. DFT+U electronic structure 

DFT+U^^'^^ is a kind of method aiming to deal with the on-site Coulomb interaction 
within the framework of DFT. It treats the local electrons in a Hartree-Fock manner, which 
drives local orbital occupations to integral occupancy as U increases. For transition metal 
oxide, one may find a larger split of d orbitals for larger U, as shown in Fig. EJ|U[ where partial 
DOS of cobalt 3d orbitals and oxygen 2p orbitals are plotted for FM, AFMl and AFM2 states 
respectively. The cobalt 3d partial DOS for Col and Co2 sites are separated. We do not 
present the results for spin restricted paramagnetic states in this section, because the energy 
differences between these states and corresponding magnetic states increase rapidly with U 
(0.4 eV/Co for U = 7.0 eV). The paramagnetic states are still metallic even with U as large 
as 7.0 eV. 

In Fig. m the top panel is the GGA result, and the three below panels represent three 
typical electronic structure behaviors for different U values. For small U (see panel b), the 
electronic structure is similar to the above GGA one, except for a strong trend towards gap 
opening and charge ordering. As shown in Fig. [3 the moments on the two inequivalent 
Co sites are nearly equal and also similar to the GGA values (which is the (f/ — J) 
limit) until U = Uc =3 eV. Above f/c, the unoccupied part of the states weight more 
and more for Co2 states with the increase of U . Disproportionation from Co^'^+ into s = ^ 
Co^^ and s = Co"^^ ions is nearly complete at [/ = 4 eV and is accompanied by a metal 
insulator (Mott-like) transition from conducting to insulating. This insulating phase results 
from the splitting of the ag states, which lead to an unoccupied narrow band containing 
one hole for each Co2 ions and an occupied band on the Col ion, as shown in Fig. We 
judge the formation of charge ordering from local magnetic moments instead of from charge 
population directly, because the charge difference between Col and Co2 are relatively small 
(only 0.1~0.2 electron). This is a result of the hybridization of the O 2p and Co 3d orbitals. 
In fact, oxygen may directly contribute to the charge ordering in some oxides, as suggested 
by Coey^S. We also notice that the charge ordering and gap opening with a moderate U is 
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far from obvious, since Zou et al}^ have gotten a metallic ground state with a relatively large 
Coulomb U (5 eV) but without adopting any supercell structure. Therefore the sodium ion 
ordering is critical to the insulating ground state. 

There is no clear boundary for moderate and large values of f/, but a relatively large 
U value will feature several characters as shown in Fig. 0Ji. First, the empty band 
are pushed further higher, which make it strongly mix with the upper Cg bands. Secondly, 
the large splitting of the Co 3(i orbitals even squeeze some Co2 "id states to the bottom 
of the O 2p manifold, and consequently the O 2p character at the edge of valence band 
near Fermi energy has exceeded the Co Sc? characters. Thirdly, the magnetic moments of 
cobalt are exceed 1 hb in this region (up to 1.42 jiB per Co for U = 7.0 eV), which seems 
curious at first sight. We now give two possible rationalizations for this observation. One 
interpretation comes from the population analysis process. We notice that oxygen is also 
spin polarized in our calculations for large U, which makes the net magnetic moment of the 
whole supercell is still about 1.0 /^b per Co. Therefore the nominal magnetic moment being 
larger than 1 fiB per Co may be just an artifact of PAW population analysis as a result 
of the strong hybridization between Co and O. Another possibility is that the large on-site 
Coulomb interaction may really promote a transition from low spin state to a intermediate 
spin state as demonstrated by Korotin et al."^^ for Co ions in vertex sharing CoOe octahedra. 
This type of transition is consistent with the strong mixing of empty ag and Cg bands. 

The evolution of electronic structures with U for AFM states is generally similar to that 
for FM state. The most significant difference is that the charge ordering occurs well before 
gap opening for AFMl state. As shown in Fig. [3 the Co^"*" moment grows immediately for 
AFMl state as U increases from 1 eV, but the gap opens until U = 3 eV, which is different 
from the results of Lee et al}^ for Naa;Co02 with x = 1/3 and x = 2/3. 

The energy difference between the FM and AFM phases as a function of U is plotted in 
Fig. 03. The general trend gives lower energy AFM state for small U and lower energy FM 
state for large U, except the relatively higher AFMl energy for small U, which is caused by 
the artificial charge ordering as previously discussed. For large U, the lower energy of FM 
state is a natural result of the weakening of superexchange effect, and it is also consistent 
with the prominent DOS of the majority spin component at deep energy as shown in Fig. HJi. 
The energy difference between FM and AFM2 is generally small according to the relatively 
small exchange interaction between different layers. According to our result, the magnitude 
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of on-site Coulomb interaction U in this material can be probed by experimental magnetic 
measurement. 

From Fig. EJini we can find an interesting splitting and ordering of the Cg orbitals. With 
the increase of U, the lower part of the Cg orbitals becomes mainly from Col sites and the 
upper part from Co2 sites. We notice that this kind of orbital ordering of unoccupied Cg 
orbitals will not affect the low energy properties of Nao.5Co02. Although it has been found 
in LaCoOs^i, orbital ordering of valence band orbitals is not very possible in this structure, 
since they are not doubly degenerated. 

A rehybridization process is proposed to describe the competition between the on-site 
Coulomb interaction and the e^-oxygen hybridization in doped CoO layers^. This process 
can be seen clearly from the partial DOS of O 2p orbitals as shown in Fig. EJEl At GGA 
level, hybridization of oxygen 2p states and Cg states are more significant than the t2g-oxygen 
hybridization. With the increase of U, the on-site Coulomb interaction becomes stronger, 
which can be minimized by unmixing the oxygen and Cg orbitals in order to decrease the 
occupation of the eg orbitals. In Fig. I^IHl we can clearly see this kind of unmixing by the 
decrease of oxygen 2p partial DOS in the Cg manifold. For large value of U, the hybridization 
of oxygen 2p states and t2g states is somewhat enhanced, which is not included in the modified 
Hubbard model study by Marianetti et alA 

We have studied the effect of Hubbard U on the electronic structures of Nao.5Co02 in the 
previous part of this section, but the actual value of U in this material is not determined 
yet, of which very different values are proposed in the literatureSiiiiSiSS. Based on the 
folowing two experimental observes, we estimate that the vaule of U should be moderate, 
as also suggested by Lee et ali^. First, the insulating ground state suggests that the value 
of U should not be too small. Secondly, Naa;Co02 with x away from 0.5 is always metallic 
indicates a narrow energy gap for Nao.5Co02, which prefer a not too large value of U . Since 
there was experimental reporli^ on electron mass enhancement for small x of Naj^Co02, 
people have compared the experimental electronic specific heat coefficient 7 with the band 
value to estimate the strength of correlationi2*i^. But in the Nao.5Co02 case, where Hubbard 
type correlation drives system to a Mott insulator, 7 will reduced by the correlation effects. 
In fact, the experimental 7 is only 3 mJ/mol-K^,^^ which is much smaller than our GGA 
band value (10.5 mJ/mol-K^), but very near to the DFT+U value (zero). 

At a plausible value of f/, 4 eV for example, intralayer AFM coupling is favored over 
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FM and interlayer AFM coupling (ref Fig. UJp)- We notice that the magnetic state in 
Na2:Co02 is far from being clearly determined in experiments. There are some experiments 
supporting AFM couplin g^°i^^ , while there are also some recent experiments supporting FM 
coupling22i2^iS. With this intermediate value of U, the ground state (AFMl) magnetic 
moment for Co2 is 1.0 fiB, and the magnetic energy is 0.123 ev/Co. 

IV. CONCLUSION 

Based on DFT+U method, we have studied the effects of Na ordering and on-site Coulomb 
correlation on the geometrical, electronic and magnetic properties of Nao.5Co02. Without 
including on-site Coulomb interaction U, the experimental insulating phase can not be 
produced by DFT. For small value of U, ground state is still metallic. At this time, no charge 
ordering occur for FM coupling, and symmetry enforced charge ordering for intralayer AFM 
coupling is unfavorable in energy. Increasing parameter f/ to a moderate value can open 
an energy gap and form charge ordering for both FM and AFM states. Large U value will 
largely split the Co 3d orbitals, which leads to squeezing Co 3d states to the bottom of O 
2p manifold, mixing of empty ag and Cg states, and large Co magnetic moments. 

Accompanying with the early work on Na2;Co02 (x = 1/3 and x = 2/3) our study 
will shed some insights into understanding of the complex interplay of Na content, super- 
structure, correlation effects and tendency of various type of ordering in Na2;Co02. However, 
a complete understanding on this issue, for example, why the x = 2/3 case with stronger 
commensurability effect and maybe also stronger correlation is not a charge ordered insulator 
in experiment, will need much further work. 
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FIG. 1: (Color online) The optimized crystal structure of Nao.5Co02. The designations of the Na 
and Co atoms are as indicated. 

FIG. 2: GGA paramagnetic band structure (left panel) and density of states (DOS) (right panel) 
of Nao.sCoOs- The Fermi energy is at zero. r=(0,0,0), M=(l/2,0,0), N=(l/2, 1/2,0), and 
A=(0,0,l/2). The dot line in the DOS panel represents the Co ag partial DOS. 

FIG. 3: (Color online) (a) Schematic diagram of the reciprocal lattice and the first Brillouin 
zone for the original P^s/mmc lattice (black) and the superstructure adopted in this work (red) 
respectively, (b)-(f) Top views of Fermi surface within the first Brillouin zone for band 96 to band 
100 of Nao.5Co02 respectively. 

FIG. 4: (Color online) Density of states (DOS) for FM phase. Total DOS, O 2p partial DOS, and 
3d partial DOS of Col and Co2 for (a) GGA result and U = (b) 3.0, (c) 4.0, and (d) 7.0 eV are 
presented. 

FIG. 5: (Color online) Density of states (DOS) for AFMl phase. Total DOS, O 2p partial DOS, 
and 3d partial DOS of Col and Co2 for (a) GGA result and U = (b) 3.0, (c) 4.0, and (d) 7.0 eV 

are presented. 

FIG. 6: (Color online) Density of states (DOS) for AFM2 phase. Total DOS, O 2p partial DOS, 
and 3d partial DOS of Col and Co2 for (a) GGA result and U = (b) 3.0, (c) 4.0, and (d) 7.0 eV 
are presented. 

FIG. 7: (Color online) Effect of the intra-atomic repulsion U on (a) magnetic moments, (b) energy 
gap and (c) energy difference between AFM and FM states of Nao.5Co02. The energy of FM state 
is set to zero. 
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TABLE I: Optimized structral parameters. Experimental values are presented in parentheses. 



Atom Site 


X 


y z 




Af 


0.003 (0) 


0.250 (1/4) 0.003 (0) 


Co2 


Ad 


0.500 (1/2) 


0.000 (0) 0.000 (0) 


Nal 


2b 


-0.030 (-0.018) 0.250 (1/4) 0.750 (3/4) 


Na2 


2a 


0.355 (0.352) 


0.750 (3/4) 0.750 (3/4) 


01 


4/ 


0.336 (1/3) 


0.250 (1/4) 0.0880 (0.0895) 


02 


4/ 


0.330 (1/3) 


0.750 (3/4) 0.0845 (0.0820) 


03 


8^ 


-0.166 (1/6) 


0.000 (0) 0.0892 (0.0895) 
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Figure 1, Li al. 
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Figure 2, Li al. 




Figure 3, Li et al. 
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Figure 4, Li et al. 




Figure 5, Li et al. 




Figure 6, Li et al. 




Figure 7, Li et al. 




